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ABSTRACT 

A recent revision of black hole scaling relations indicate that the local mass density in hlack holes should be five times higher than 
previously determined values. The local black hole mass density is connected to the mean radiative efficiency of accretion through 
the time integral of the AGN volume density and a significant increase in the local hlack hole mass density would have interesting 
consequences on AGN accretion properties and demography. One possibility to explain a large hlack hole mass density is that most 
of the hlack hole growth is via radiatively inefficient channels such as super Eddington accretion; however, this solution is not unique. 
Here we show how it is possible to accommodate a larger fraction of heavily buried, Compton-thick AGNs, without violating the limit 
imposed by the spectral energy density of the hard X-ray and mid-infrared backgrounds. 
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1. Introduction 


P 

strong correlations between the masses of the central super 
0 massive black holes (SMBH) and the global properties of their 
^ host spheroids such as luminosities, dynamical masses and ve- 
c/1 locity dispersions, i.e. the scaling relations, can be used to con- 
, ^, vert the mass function of local galaxies into a black hole mass 
function and, by integration, into a local SMBH mass density 

P»Joc- 

A recent comprehensive analysis of black hole mass mea- 
surements and scaling relations concluded that the canonical 


(N 


black-hole-to-bulge mass ratio, MsMBHlMsuige, shows a mass 
dependence and varies fro m 0.1-0.2% at Mbuige - 10 ^Mq to 
~ 0.5% at Mbuige = 10'*Mo (Graham & Scott ^20 13 Kormendy 


j |& Ho |[2013| l. The revised normalization is a" factor of 2 to 5 
larger than previous estimates ranging fro m ^ 0.10% ([Merrit & 


|Ferrarese [2001 McLure & Dunlop [2002 Sani et al. 20111 to 


X 


2008 |l. Assuming a constant radiative efficiency and neglecting 


mergers, the continuity equation can be integrated to obtain 


.C" = Ur 


l-e 


- {kBOlWx 


l-e 


( 1 ) 


where e is the radiative efficiency, Ut is the bolometric, comov¬ 
ing energy density of AGNs described by a bolometric luminos¬ 
ity function <1)(L, z). 


Ur — 


MJ 


L<1>(L, z)dL. 


( 2 ) 


0.23% ( Marconi & Hunt ||2003| l, therefore resulting in an over¬ 

all increase in the normalization, which is dominated by massive 
bulges. 

The mass function expected from black holes grown by ac¬ 
cretion during active phases and the corresponding integral over 
the cosmic time p.^acc, are closely related to the AGN luminosity 
function via a mass-to-luminosity conversion factor. The argu¬ 
ment, originally proposed by |Soltarr| ( |1982[ ), can be formulated 
as a continuity equation for the mass density of black holes (e.g. 
Marconi et al. |[2004l IShankar et al. |[2004| [Merloni & Heinz | 


The comoving energy density Ut can be written as the prod¬ 
uct of the X-ray energy density Ux times an average bolometric 
correction factor {ksoL) to convert the X-ray emissivity to bolo¬ 
metric emissivity (e.g. jLusso et al. |2012| l. 

Equation 1 is widely used in AGN demographic studies to in¬ 
fer the efficiency of accretion processes by comparing the values 
obtained by the integration of the bolometric luminosity function 
with the local black hole mass density obtained from the scaling 
relations. A good match between p.^ioc and p.^acc is achieved for 
an accretion efficiency consistent with that expected from a stan¬ 
dard Shakura & Sunyaev ( |1973 1 accretion disk, broadly sup¬ 
porting the possibility that most of the SMBH mass is efficiently 
accreted over the Hubble time (e.g. [Eabian & Iwasawa~||1999| l. 
Using the X-ray background as an integral constraint and as¬ 
suming that the bulk of the X-ray emissivity is at (z) ^ 2, [Elvis 


[et al. [ ( 2002| l infer a lower limit on the accretion efficiency of 
e >0.15, arguing that most SMBH are rapidly spinning. 

Most recent estimates are in good agreement with (e) 0.1. 
Eor instance, Marconi et al. [ ([2004)1 use d luminosity functions 
in the optical 3-band ([Boyle et al. [2000 1, soft X-ray, 0.5-2 keV 
band (Miyaji et al. 2000| l, and hard X-ray, 2-10 keV band ( [Ueda[ 


[et al. 


2003|l transformed into a bolometric luminosity function 


using a luminosity dependent X-ray to optical spectral index 
aox ( [Vignali et al. [|^03[ ). The correction for the number of 
missed AGNs in hard X-ray surveys (i.e. Compton-thick) is ob¬ 
tained by requiring that the observed X-ray luminosity function 
matches the X-ray background spectrum. The resulting estimate 
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of p,_acc, assuming e=0.1, is fully consistent with the lo cal mass 


density p.joc- 4.6 x 10 ^Mq Mpc More recently, 


-3 


Ueda et 


al. (2014[l published an updated version of the hard X-ray lumi- 


nosity function, including an estimate of the Compton-thick pop¬ 
ulation needed to fit the hard X-ray background spectrum. They 
reached similar conclusions on the average accretion efficiency 
Ho^ins_et al. | ( 2007[ l bolometric corrections and the 
2009 1 best fit value for p.joc- 4.9+0.7 x IO^Mq 


using the 


Vika et a 
Mpc^. 

The agreement of the two approaches described above is not 
surprising despite the different assumptions of the luminosity 
functions and bolometric corrections. In fact, both are linked to 
reproduce the hard X-ray background, which can be considered 
an integral constraint on the total mass accreted over the cosmic 
time and locked in SMBH. The almost identical values of the as¬ 
sumed local density {p.joc) thus lead to very similar conclusions 
on the average accretion efficiency. 

At face value, and keeping all the factors entering in Eq. 1 
fixed at the fiducial values described above, a higher value of 
the integral local mass density would lead to a similarly lower 
value for the average radiative efficiency. The possibility that ra- 
diatively inefficient accretion processes may be more common 
than previously thought has been put forward by |Novair| ( |2013| l. 
His argument is based on a simple linear scaling of the com¬ 
monly accepted value for the accretion efficiency (e ~0.1) by a 
factor of 5, the maximum discrepancy between the revised esti¬ 
mate ( [Graham & Scott [2013 [Kormendy & Ho [2013| l and previ¬ 
ous values. Radiatively inefficient processes (i.e. slim accretion 
disks) are advocated to explain the fast growth of SMBH in the 
early Universe (e.g. jMadau et al. |2014| l. 

The revised estimate of the local mass density of SMBH may 
have important consequences for AGN demography and accre¬ 
tion physics. In the following we discuss the impact of an in¬ 
creased fraction of heavily obscured SMBH with respect to cur¬ 
rent observational constraints. 


2. A large population of deeply buried SMBH? 

The contribution to the SMBH mass density of unobscured and 
obscured AGNs and the degeneracies intrinsic to the comparison 
between p.^acc and p.joc, are manifest if equation 1 is written as 


2 1 — e 

P.,accC - {kBOl) -Uxo(l + / J RqBs), 


(3) 


where Uxo is the X-ray energy density of unobscured AGNs, 
as determined from the soft X-ray luminosity function (e.g. 
Hasinger et al. 2005| l, while {kboi) is the average bolometric cor¬ 
rection which is assumed to be the same for obscured and unob¬ 
scured AGNs (e.g. Lusso et al. |2012| l. The sum is over the ratios 
to unobscured AGNs of the obscured population contributing to 
the SMBH mass density and the X-ray background; 


OBS - RThin+RThick - Rrhin + RMThick + RuThickX fs). 

(4) 


In the following we refer to the Marconi et al. j ( |2006| ) cal¬ 
culations, which are based on the Gilli, Comastri & Hasinger 
(2007, GCH07) AGN synthesis model where the population of 
Compton-thick (CT) AGNs is split into mildly CT (MCT) when 
the optical depth for Compton scattering tj is of the order of 
a few, and heavily CT (HCT), or reflection-dominated, when 
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tt" » 1. The former provide the largest contribution to the XRB 
(X-ray background) spectrum at 20-30 keV, while the contri¬ 
bution of the latter depends mainly on the fraction of the repro¬ 
cessed radiation fs reflected or scattered towards the observer 
which, in turn, depends on the geometry of the absorber. In 
GCH07 the fs value in the 2-10 keV band is assumed to be 
0.02, in broad agreement with the observations of Compton- 
thick AGNs in the local Universe (see Fig. 1 in GCH07). More¬ 
over, the ratio between Compton-thin and unobscured sources 
is luminosity dependent and smoothly decreases from R = 4 at 
low luminosities (logLx = 42) to R = 1 at high luminosities 
(logLx = 46). The luminosity dependent obscured fraction is 
assumed to be the same for both Compton-thin and Compton- 
thick objects. As a result, the total (thin plus thick) obscured- 
to-unobscured ratio decreases from 8;1 at low luminosities to 
2:1 at high luminosities. The luminosity averaged ratio is ap¬ 
proximately 3 which means that for each unobscured AGN there 
are three Compton-thin and three Compton-thick sources (1.5 
MCT and 1.5 HCT) and thus 'ZRobs—(>- An excellent match be¬ 
tween p. acc and p. loc, satisfying the XRB constraints, is found 
for e =0.i ( [Marconi et al. 12006) 1. 

We consider the most conservative case of a local mass den¬ 
sity which is a factor of 2 higher than previously considered. The 
revised value ( [Graham & Scott [2013[ Kormendy & Ho [2013 1 is 
almost exactly a factor of 2 higher than that obtained by Marconi) 
[et al. [ ( [^04] l and [Vika et al. [ ( [2009l l. We also assume that the ad¬ 
ditional mass density is not due to a lower efficiency, but to a new 
population described by the same global parameters (cosmolog¬ 
ical evolution, bolometric correction, and radiative efficiency) of 
the known AGNs responsible for the XRB. If this were the case, 
equation 3 becomes 


2p 

• ,acc^ - {kBOL}UxO -(1 + 'y'^RpBS + Rnew), (5) 

/ 

where Rmew is the contribution of the new population respon¬ 
sible for the local mass density excess. It is straightforward to 
show that 

Rnew - {^ + ^^Robs) - (6) 

implying that the contribution to the mass density of the new 
population is the same as that already contributing to the XRB. 
It goes without saying that the new population would vastly ex¬ 
ceed the limits imposed by the XRB itself, unless the sources 
are so extremely obscured that X-rays are almost completely 
suppressed. In other words, the upward revised normalization 
can be matched if each SMBH contributing to the XRB, has an 
X-ray quiet counterpart that adds mass, but almost no X-ray 
radiation. 

This figure could be mitigated noting that the present uncer¬ 
tainties on the XRB spectral intensity at ^ 20-30 keV are such 
that it would be possible to accommodate additional sources with 
the same spectrum of HCT in GCH07. More specifically, if the 
reflection-dominated sources are four times more numerous than 
in GCH07 (Fig. 1), the predicted XRB flux at 20-30 keV would 
still be consistent with the data, within the measurement errors. 
If this were the case the size of the X-ray quiet population would 
be reduced to R^^j^ = 2.5. This new value is derived by the 
equality 


Rhct + Rnew - 4 x Rhct + R 


NEW’ 


(7) 
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Fig. 1. The broadband cosmic XRB spectrum. The thick, solid, dark blue curve and the thin, solid, black curve show the total AGN spectrum and 
the contribution of Compton-thick (mildly-theavily) AGNs in the GCH07 model, respectively. The thick, dashed, dark blue curve and the thin, 
dashed, black curve are as above, but assume a 4 times larger abundance of heavily Compton-thick AGNs. The XRB measurements from various 
mission/instrument combinations are marked by different symbols as labeled. Labels on the left (right) generally refer to measurements at £ < 10 
(> 10) keV. The relevant references are as follows. ASCA GIS: Kushino et al. (2002); ASCA SIS: Gendreau et al. (1995); ROSAT: Georgantopoulos 
et al. (1996); BeppoSAX/MECS: Vecchi et al. (1999); XMM(2002): Lumb et al. (2002); XMM(2004): De Luca et al. (2004); RXTE: Revnitsev et 
al. (2003); Swift/XRT: Moretti et al. (2009); Nagoya balloon: Fukada et al. (1975); SMM: Watanabe et al. (1997); HEAO-l: Gruber et al. (1999; 
open circles and filled triangles), Kinzer et al. (1997; hatched area); BeppoSAX/PDS: Frontera et al. (2007); INTEGRAL(2007): Churazov et al. 
(2007); INTEGRAL(2010): Turler et al. (2010); Swift/Bat: Ajello et al. (2008). 


we note that Rhct is 1.5. In this scenario, only about 18% of 
the SMBH population would be X-ray quiet, but having sig¬ 
nificantly increased the fraction of HCT sources in the GCH07 
model, the fraction of deeply buried AGN {HCT plus X-ray 
quiet) would be about 60% of the total SMBH population. 

We can relax the GCH07 assumption of the reflection yield 
{fs - 0.02), and assume that both the GCH07 model HCT 
sources and the new population are described by a heavily ob¬ 
scured, reflection-dominated spectrum, with the same value of 
the reflection yield in the 2-10 keV band. In this case the new 
values of the yield of reflection-dominated sources are obtained 
from 


Rhct R 0.02 - f^Ew ^ [^hct + ( 8 ) 

where Rhct-^-^ is the fraction assumed in GCH07. With 
these assumptions, we obtain /Ar£w=0.0035 and for 

Rnew-'^ or Rnew -'^-^7 depending on the assumed level of the 
XRB. 


3. Discussion 

A large population of heavily obscured, reflection-dominated 
AGNs is considered here to All the mass gap between the most 
recent estimate of p.joc and previous values. In order to accom¬ 
modate a significantly higher number of SMBH without violat¬ 
ing the limit imposed by the X-ray background, the hypothetical 


new population should be characterized by an extremely faint, 
almost vanishing, X-ray emission. 

The suppression of the hard X-ray emission (up to a few 
hundred keV) would require both extremely high column den¬ 
sities (t » 10) and large covering factors (approaching 47r). 
The detection and the measure of column densities in the high 
T regime {Nh ~ 10^^ cm“^ or even more) is extremely challeng¬ 
ing. Large covering factors of the obscuring gas, with reflection 
yields much lower than 1%, are invoked to explain the obser¬ 
vations of a few heavily obscured and Compton-thick AGNs 
observed with Suzaku (e.g. Ueda et al. ||2007[ Comastri et al. 


2010|l. There is also tentative evidence for a larger covering fac- 


tor for increasing absorption among Chandra Deep Field South 
faint sources ( |Brightman et al. ||2014| l. The upper limits on the 
2-10 keV X-ray emission of a small sample of luminous AGNs, 
hosted by nearby ULIRGS ( [Nardini & Risaliti~||2011| l, are so 
tight that the corresponding upper limits on /j are of the order of 
10“^ - 10 If the presence of these AGNs is widespread among 
infrared galaxies, they could be excellent candidates for the new 
population. 

At face value, a population of reflection-dominated 
Compton-thick AGNs, with a suitable combination of covering 
factor and column density, would still be consistent with the up¬ 
per envelope of present uncertainties of the hard XRB spectrum. 

While the new population could have so far escaped direct 
detection in present deep, soft (< 10 keV) X-ray surveys, and 
would remain largely undetected in shallow NuSTAR hard X- 
ray surveys ( [Alexander et al. ||2013| l, the reprocessed AGN lu- 
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minosity will emerge in the infrared. A first order, integral esti¬ 
mate may be obtained using the hard XRB spectral intensity as 
a limit of the total energy output combined with suitable bolo- 
metric corrections (e.g. [Fabian & Iwasawa~||1999| [Risaliti et al.j 
[] |2002| ). The results indicate that the AGNs responsible for the 
XRB make 10-15% of the IR background around its ~100/rm 
peak. Their contribution increases at shorter wavelengths, be¬ 
ing of the order of 20^0% at 60/rm and possibly higher at 


15/rm depending on the dust effective temperature (Risaliti et 


jal. [2002] ) and the relative fraction of Compton-thick sources in 
the adopted XRB model. The impact of the new population of 
heavily Compton-thick AGNs hypothesized here is even more 
uncertain. In the most naive approximation, the above estimates 
should be rescaled by a factor of 2 making the AGN contribution 
dominant at short wavelengths (< 20 pm). However, according 
to a recent, comprehensive population synthesis model of the X- 
ray and infrared backgrounds ( |Shi et al. |2013| ), the contribution 
of the X-ray emitting AGNs to the 100/tm background is negli¬ 
gible and is of the order of 15% at 5-10/im. If this were the case, 
there would be plenty of room to fit a large population of heavily 
obscured SMBH without exceeding the limits imposed by the IR 
background. 

Observational constraints on the abundance of sources likely 
to harbor a deeply buried new AGN may be obtained by mid- 
infrared photometric and spectroscopic surveys carried out with 
the Spitzer Space Telescope. Highly obscured and Compton- 
thick AGNs are likely to be affected by heavy dust absorption 
and steep mid-infrared to optical slopes. A widely adopted se 


lection method is based on the optical to lApm color (e.g. Fiore 
jet al. ||7009| ). The optical depth of the silicate absorption fea¬ 
ture at 9.1pm is also considered a reliable tracer of heavy dust 
absorption. Good quality, mid-infrared spectra were obtained for 
local luminous infrared galaxies in the GOALS survey (jStierwaltj 
jetal. [201 3| l and z -0.7 COSMOS galaxies ( |Fu et al. [201 0| l. The 
fraction of optically thick (T9.7 >1) galaxies, which are likely to 
host a buried AGN on the basis of SED decomposition, is es¬ 
timated to be of the order of 10% or even lower. However, it 
should be noted that current mid-infrared photometric and spec¬ 
troscopic surveys are relatively shallow and, moreover, there are 
many Compton-thick AGNs that do not display strong silicate 


absorption features (e.g. Georgantopoulos et al. 2011 1 . Given 


the present uncertainties and the difficulty of distinguishing the 
hypothetical population of deeply buried sources from normal 
and star-forming galaxies, we conclude that the new population, 
is unlikely to exceed the IR background spectral intensity. Can¬ 
didates of the new population may be found among sources with 
a steeply rising SED from the near- to the far-infrared and/or 
among those with the highest extinction in the silicate absorp¬ 
tion feature at 9.1pm. 


4. Conclusions 

We put forward the possibility that a sizeable population of hith¬ 
erto uncovered ultra obscured AGNs would explain the revised 
upward value of the SMBH local mass density p,joc- 

We are aware of the strong degeneracies between the var¬ 
ious parameters used to build our argument and, in particular, 
the radiative efficiency e. Within the reasonable assumptions that 
the evolution of the X-ray luminosity function and bolometric 
corrections are robustly constrained by recent surveys, and the 
spectral intensity of the X-ray background and the accretion pro¬ 
cesses are self similar and not evolving, a large number of deeply 
buried SMBH may be a viable solution beyond, radiatively inef¬ 
ficient accretion, for example. 
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This would imply that the column density distribution of nu¬ 
clear obscuring gas in SMBH is top heavy and skewed towards 
values of the order of cm“^ or even higher and the cover¬ 

ing fraction of obscuring material may approach Ttt. Binding ob¬ 
servational proof of the existence of these ultra-obscured SMBH 
is challenging. Deep mid-far-infrared photometric and spectro¬ 
scopic surveys, matched with sensitive X-ray surveys would be 
needed to estimate the fraction of infrared emitting AGNs with¬ 
out an X-ray counterpart and thus presumably heavily buried. 
High density gas illuminated by a hidden nuclear source may 
be traced by molecular transitions observable in the far-IR with 
ALMA ( jlmanishi & Nakanishi |2014| . A fraction of them could 
be found among the sources responsible for the unresolved XRB 

2012| l. We are hopeful that a 


in the 6-8 keV band (Xue et al. 


few examples will be uncovered by forthcoming deep hard X- 
ray surveys with NuSTAR ( [Harrison et al. |2013| l and ultra-deep 
ATHENA surveys. 
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